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ABSTRACT 


Multlspectral analysis of LANDSAT imagery provides a rapid 
and accurate means of identification, classillcation, and meas- 
urement of strip-mined surfaces in Western Maryland. Four 
band analysis allows distinction of a variety of strip-mine associ- 
ated classes, but has limited extendibility. A method for surfaee 
area measurement of strip mines, which is both geographically 
and temporally extendible, has been developed using band-ratioed 
LANDSAT reflectance data. The accuracy of area measurement 
by this method, averaged over three L \NDSAT scenes taken 
between September 1972 and July 1974, is greater than 93%. 

Total affected acreage of large (50 hectare/120 acre) mines can 
be measured to within 1.0%. 



LANDSAT INVENTOI Y OF SURFACE- MINED AREAS 
USING EXTENDIBLE DIGITAL TECHNIQUES 


I NTRODl'CTlO F 

Extraction of ncar-surface coal by surface or "strip" mining techniques is an 
economically feasible response to the Nation's growb demand for enei^. In 
order to reduce dependence on foreign imports, co»I mining activity has increased 
greatly in the last year, and will continue to do so as new reserves are opened. 
The Nation's need for energy, however, must be balancc-c with the environmental 
consequences of coal exploitation. Surface mining unaccompanied by reclamation 
renders the land useless for other productive uses. Property near or adjoining 
mine sites is degraded in value (Ref. 1), and severe erosion, landslides, flood- 
ing, ait and water pollution may also occur. 

Of the estimated 2. 2 million acres of land stripped for coal in the United States, 
only a third has been reclaimed (Ref. 2). Frequent monitoring by regulatory 
agencies is required to insure reclamation success. National reclamation stand- 
ards do not yet exist and requirements vary widely among the 23 states which 
presently regulate strip mining. Thus, information on the location, size, and 
condition of mines is often lacking or inadequate. In a cooperative NASA/State 
of Maryland effort. Imagery relayed from the NASA Earth Resources Technology 
Satellite (LANDSAT-1)' was applied to the monitoring Information needs of the 
Maryland St^te Bureau of Mines. The objectives of this cooperative study were 
(1) determine the accuracy of satellite data for measuring strip mines of the 
size common in Western Maryland, fjid (2) develop an operationally feasible 
procedure for large area inventorying and monitoring of surface mining. 

BACKGROUND 

The State of Maryland became officially involved in the regulation, monitoring 
and reclamation of strip mines in 1967 when the State's coal strip mining law 
w’as enacted. This law, amended in 1969 and 1971, established a Land Reclama- 
tion Committee which administrates the provisions of the act. Under this act, 
strip mine operators are required to obtain licenses and permits, submit mining 
and reclamation plans, procure performance bonds, and periodically report on 
the amount of land area disturbed. All of the disturbed surface, including storage 
areas for topsoil and spoils, haul roads, and areas disturbed by the movement 
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of equipment, must bo reclaimcHi aecordlnK to State -approved olans. Mines 
closed before the law's enactment will be reclaimed by the Stat J. 

To assist in the tasks of planning, monitoring, enforcement, and reclamation, 
the Maryland State Ikireau of Mines needs information on the location, size, 
and condition of coal surface mining areas. Three mine inspectors presently 
report monthly on active mining operations within a 15-11 square kilometer area 
(595 square miles). The new regulations and increased coal production \\111 
require additional inspectors and more efficient means of data collection. 

STUDY ARLA 

The region studied in this investigation includes most of Garrett County and the 
adjacent area of Allegany County in Western Maryland. The State's only coal 
deposits are located in these counties. This region, part of the Allegheny 
Plrteau of the Appalachian Itange, consists of gently folded mountains, the result 
of differential erosion of sedimentary materials which were originally deposited 
in horizontal layers and were compressed to form a series of parallel mountains 
and synclinal basins. The coal-bearing strata are of Pennsylvanian age and lie 
within 200 feet of the surface in five major basins (Figure 1). The easternmost 
basin, Georges Creek, contains 233 square kilometers (90 square miles) of 
intensive strip mine activity. Most of the strip mines in this area follow topo- 
graphic contours and are 75 to 400 meters wide (250 to 1300 feet). High walls 
range from 20 to 45 meters (60 to 150 feet). The majority of mines are less 
than 90 hectares (220 acres) in size. It was within the Georges Creek basin 
that multispectral training of LANDSAT data was conducted. 

METHODS 

AREA MEASUR^.MENTS 

Accurate area measurements from aerial photographs were required for this 
study in order to quantitatively inventory the r.rea affected by surface mining 
and verify the accuracy of the satellite Inventory. Prior to 1975 the State 
recorded only the area of coal annually exp>osed. For this study, there- 

fore, total affected area was determined by planin '^try of NASA/Wallops low 
altitude aerial photography taken in October 1973 (Figure 2). The area measure- 
ments obtained from these photographs were used to verify the accuracy of the 
LANDSAT classifications for eight surface mines in the test site. 

Obtaining accurate area measurements from low altitude aircraft Imagery is a 
difficult, time consuming task due to a combination of factors: off-nadir viewing. 
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the varying elevation of the mines in this area, and the distortion due to the 
mountainous terrain. The following procedure was followed in calculating the 
area of each test mine from a.'rcrnft photography: 

• The mine areas were planimetered four times with precision among 
these being greater than 96 %, (Readings with a large deviation were 
rejected. ) The average of these readings was used in the calculation. 

• A table of photo scale versus aircraft altlbide was calculated using 
the formula: 

Photo scale » 

• The table of photo scale was verified and corrected using ground 
distance measurements on the aerial photographs and topographic 
maps. 

• The camera height was determined by subtracting the mean topographic 
elevation of the mines from the corrected aircraft altimeter readings. 
Large mines were segmented into sections of nearly equal elevation. 

V Using the appropriate photo scale factor the planimetry values were 
converted to acres. 

• The mine acreage values were corrected for the viewing aspect angle. 

The surface area measurements obtained in this maiaicr are shown in Table II 
which is in the Results Section of this paper. The test mines shown in Figure 2 
range in size from 12.7ha (31.4 acres) to 98.7ha (243.8 acres). Surface area 
values from aircraft photos are 25 to 30 percent larger than the publisheJ • 
figures wh’ch at that time related only to the coal surface exposed; the ail. . • 
values are for the total affected area including spoil piles and haul roads, ;* 
accurate reclamation projections, the State now requires data on this total dis- 
turbed area. 

High altitude aerial photographs of the Georges Creek basin are available from 
NASA/Ames for 1972 rind 1974. One of the test mines, Franklin Hill-A, was 
planimetered on these photogra()hs because the entire mine was not included on 
the 1972 low altitude photos. In general, however, the areas of interest were 
too small and too near the format edge to permit accurate planimeter measure- 
ment. 


cal length 
Altitude 


3 



The high altitude photography did show that most of the test mines changed very 
little in size during the two-year period between the photos. The only tet i mines 
whose areas were significantly different from year to year arc those located 
near Mill Run. Mill Run-B had not been opene<l in 1972, but was completely 
stripped by October 1973. In 1974 it was backfilled. The lower section of 
Mill Run- A was bare in 1972 and almost completely revegetated in 1973 (sec 
Figure 2). The other test mines remained approximately the same size from 
1972 through 1974. 

MULTISPECTRAL ANALYSIS 

An interactive, multlspectral image analysis system^ was used to perform the 
digital data processing of LANDSAT computer compatible tapes (CCTs). The 
system's image analysis console houses a color image display, controls, and 
a special purpose high speed processing logic. A mini-computer together >\1th 
peripherals serves as a system process ccmtroller anu computational device. 
The user interacts directly with the computer through a graphics entry/display 
terminal. The graphics terminal also serves to display the quantitative proc- 
essing results in both numerical and graphic form. 

Four- Band Classification 


A portion of LANDSAT-1 scene 1405-15242 (1 Sopt. 73) containinr the Georges 
Creek basin was digitally enlarged to fill the system's 512 X 512 picture element 
color display (Figure 3). The area displayed is approximately 51.8 square 
kilometers (20 square miles). Supervised training and classification were then 
performed on the four LANDSAT spectral bands. Training sites were selected 
with an electronic cursor which is sized and positioned using a Joy-stick. The 
image analysis system's special purpose hardware Identifies the spectral reflec- 
tance range within the training site in the four LANDSAT bands simultaneously. 
The minimum and maximum reflectance values in each channel (band) of the 
training area are then used to define the limits of a 4-dimensional spectral paral- 
lelepiped. The picture elements of the entire displayed image are examined 
pixel-by-pixel. Those pixels lying spectrally within the parallelepiped bounds 
defined by the training site are identified or "alarmed" on the TV monitor. The 
entire process requires less than 5 seconds. The system user then has the 
option of modifying the spectril signature to increase or decrease the alarm 
through thresholding the parallelepiped boundaries. This procedure is known as 
single parallelepiped training and classificatio’'- 


^ r^ral Electric IMAGE 100 System 
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In this study s.'ngle-paral'elcpipcd training was applied to areas within the 
Georges Creek basin whose surface cover was known from interpretation of 
the low altitude aerial photography and frr-,n field inspection. In certain eases 
pixels were classified into more than one category. Aerial photographs were 
then consulted to Hclcrmine the proper classification, and the signatures of the 
overlapping clashes were modified to eliminate the conflicts. 

Seven classes, five of w'hich describe areas affected by strip mining, were 
identified. Discussion of these classes and their extension throughout the study 
region is presented in the Results Section of this paper. 

itand- Ratio Classification 


Tkie great variability among 4-band signatures for various types of strip-miniiig 
su /faces discourages the use of this means of classification on an operational 
monitoring basis. A great deal of ground verification data is required, and train- 
ing must be carried out on each type of surface. This is not ur'v true in Western 
Maryland, but the Environmental Protection Agency (Ref. 3) has reported great 
variation among signatures for strio mines studied in W'yoming and Montana. 
Preprocessing of the LANDSAT multispectral data before classification provided 
a means of distinguishing strip mines in a single classification. The objective 
of preprocessing is to transform the sensor outputs to minimize the effects of 
environmental, observational and sensor conditions on signature extraction, 
exher investigators (Ref. 4) have found that band-ratioing techniques minimize 
systemic errors and decrease temporal and geographical differences. The 
rationale behind ratioing can be illustrated by a simplified model of spectra' 
signal in the narrow bandwidth (i) of a sensor: 


Sj' * m, Sj + aj 


where Sj' is the observed reflectance and Sj is the reflectance at the surface. 
If the multiplicative terms are larger than the additive terms and these are 
approximately invariant over adjacent spectral bands, then when i+1 = J, the 
ratio 


£l . 

m,^ = 

Sj 

Sj 

mjSj 

Sj 
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Is not affected by the multiplicative error factors. Hence, multiplicative system 
errors are minimized In the ratioed sl|{nals. Table I, borrowed f*X)m Krlegler 
(Itef. 5), Indicates the factors that produce si(^al vanaticn. Note that most of 
the factors related to the useful si^^nal arc multiplicative. On the other hand, 
the effects of atmospheric backscatter and noise arc additive, and they tend to 
mask the useful si|{nal. 


Table 1 

General i'acturs That Produce Signal Variations 



Causes and Dependencies 

Type of Factor 



MulU. 

Add. 

Illumination 

Shadows, time of day, clouds, etc. 

X 


Transmittance 

Altitude, liaze, aerosols, scan angle 

X 


Ueflectance 

Scar angle, sun angle, species, 
maturity, vigor 

X 


Atmospheric 

Backscatter 

Altitude, ha/e, aerosols, sun angles 


X 

Sensor Gain 

Different setups, different days 

X 


Noise 

System components 


X 

(primarily) 


Identical portions of the 1 Sept. 73 (1405-15242) and the 6 Sept. 72 (1045-15245) 
images were prepr^ cessed on the multispec.tral image analysis system to yield 
twelve band- ratios for each date. The procedure used to determine the most 
accurate strip mine classification using these ratios was as follows: 

• In both Images polygon-shaped training sites were defined to surround 
each test strip mine. 

• Histogram lists of the polygon areas were obtained for all twelve 
LANDSAT ratios. These provided pixel counts for each mine in 
several ^*ray level ranges. 
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• Pixel counts for each mine were conve’’ted to surface area. A LANOSAT 
pixel equals 0.453 hectares (1.12 acres). 

• A least squares regression was then ralculated for the aircraft area 
values (y) and the LANDSAT area vr-’. es (x) for each gray level range, 
using data from both 1972 and 1973 imagery. The regrcs«ion result- 
ing in the smallest standard error of estimate was selected as provid- 
ing the most consistently accurate surface area measurement from 
both 1972 and 1973 data. 

The smallest standard error of estimate, 2.16 hectares, was obtained using 
the bandrat'.j of MSS 5/6 (Figure 4). The resulting equation was: 


y ■ 0.295 + 0.970 x + e 


where y is the "true" or aircraft area in hectares; x is the LANDSAT area in 

hectares; and ( is the residual. This equation was used to "adjust" the LANDSAT 

area measurements for the two images. The MSS 5/6 strip mine signature 

selected by Inis procedure was applied to a third satellite image, 1729-15164 of 

22 July 1974. The pixel counts for each test mine were converted to hectares 

and adjusted by the factors in the regression equation. The resulting surface I 

area values for all three ima,;es are given in Table II. 

RESULTS 

FOUR BAND CLASSIHCATION 

Using the multispectral analysis system, 4-channcI training and classification 
on sites of known surface characteristics resulted in the identification of seven 
classes; three strip mine surface classes, bare soil, partially revegetated, and 
two classes of undisturbed land including forest and open fields. Figure 5 
demonstrates the results of this classification of the test site on 1 September 
1973 LANDSAT-1 imagery. 

Field analysis and communication wltli the Bureau of Mines revealed several 
additional details about these classes. 

• Strip Mine Class 1 corresponds to the exposed subsoil and the spoil 
piles of relatively light color, particularly from the Franklin and 
Barton coal seams. The mines in this class were either open or 
backfilled with spoil material at the time of the satellite overpass. 
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• Strtp Mine Class 2 det.'^es the open or backfilled mines of darker 
subsoil, specifically ti.ose associated with the Pittabur^t coal seam. 

The ovcrbui Jon of this seam contains a large quantity of burned shale 
(rashlngs), which has a characteristic red color, and therefore, a 
distinct reflectance signature. The red shale Is often used as road 
surfacing material by the mine operators, and some roads are Included 
In this class. 

a Strip Mine Class 3 corresponds to strip mines that have been backfilled 
and graded, some of which are being used os landfill areas. Most of 
these nilnes cxpoeetl the L'ppcr Hakerstown coal seam and Its associ- 
ated spoils. 

a The Bare Soil Class Identifies strip mine surfaces which have been 
backfilled with the spoil, graded and covered with topsoil. The soil 
Is brighter than subsoil material In all four LANDSAT spectral bands 
and therefore, this class Is distinct from the three described above. 

Some of the mine surfaces Included In this class have been seeded but 

are not yet revegetated. • 

a The Revegetated Class relates to surfaces which arc similar In cover 
to the sparsely grassed-over airstrip on Franklin Hill, which Is an 
old strip mine that has been approved as revegetated by the State. The 
boundaries between open strip mines and forests or fields are some- 
times falsely Included In the revegetated class. This Is due to the 
fact that the average reflectance of the flcid of view of these boundary 
pixels Is similar to that of snarsc vegetation on closed strip mines. 

a Open Field refers to unforcsted land, chiefly crop and grazing land, 
which has never been s'.iipped. 

• Forest refers to areas covered by dense trees and scrub vegetation. - 

The signatures derived for these seven classes were applied to the entire study ' 

area which contains the adjoining areas of Garrett and Allegany Counties In 
Maryland, and Grant and Mineral Counties In West Virginia. The resulting 
regional classification Is shown In F'igure 6. 

Certain surfaces are falsely alarmed (classified) when the multlspectral signa- 
tures are extended over a larger area. For example, the town of Keyser, 

Maryland, located to the right of center In Figure 6, Is classified as strip mine. 

Apparently the vegetated/non-vegetated contrast In reflectance of the town Is 

similar to that of the mines. Another false alarm Includes the railroad yard 

west of Keyser. The dark railroad bed material Is classified in the second I 

strip mine class. The bare soil classification correctly Identified areas other 
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than backfilled spoils. Alung the west bonk of the North Hraneh of the IHitomac 
lUver, tlie construction site of a new rail cut lot'uted adjacent to the river was 
classified as bare soil. 

BANIVHATIO C1.ASSIFICATI()N 

rsinK band-ratioed data and linear regression :uiulysis it was possible to extend 
a single signatun* to three LANDSAT images covering a time span of two years. 

Of the 12 LANDSAT ban<l- ratios, the ratio of MSS 5/6 proved to provide the 
most consistently accurate results when comparcxl to aircraft plonimi'try. Using 
the linear regression equation the standard error of estimate is less than 5 
pixel-sl/ed units. The correlation coefficient (r) of the linear regression is 
0.997. The results of the band-ratio classification of the Sept. 72, Sept. 73, 
and July 74 LANDSAT data arc presented in Tal)lc li and Mgure 7 and 8. 

The average difference from the aircraft photography is ±6.9',( for all three 
years. The largest errors occur in the July L974 image; six of the eight mines 
arc underestimated by an average of 121.. Comparison of the 1974 classification 
with the 1973 and 1972 classifications and with available air photos suggests that 
most of this ch •*<;, •* is due to natural revegetation around the edges of inactive 
mines. Tht; ' n •• .ocnix liill mines (I-'igure 2) ceased operations l)eforc the 
passage of uie Maryland Strip Mining Bill in 1967. They are void of topsoil, 

.vid natural revegetation is * regressing slowly. This natural recovery is seen 
in the coctinued decrease c the 1972 thro igh 1974 satellite area measurements. 
The two Aaron Run mines arc also revegetating. These mines were forfeited 
to the State unreclaimed prior to 1972. The State has backfilled and planted 
these areas, but the LANDSAT data indicates that reclamation has not been as 
burcossful us in other parts of the test site. 

:cate records indicate 137 acres of coal were exposed and 137,000 tons of coal 
removed from the Mill Run-A mine (Figure 2> during the perio<l monitorexl by 
the LANDSAT images. Aerial photography taken in October 1973 shows that 
the total affected area was 217.9 acres, 14. 5 acres of which was partially 
reclaimed. The Sept. 73 LANDSAT classification identified 198.3 acres of 
unreclaimed area affected by surface mining, a 2. 6l difference from the aerial 
interpretation. More reclamation had occurred by the time of the 1974 LANDSAT 
image which shows 191.7 acres in Mill Run-A. The 1972 LANDSAT image shows 
the full extent of the mine (222.1 acres) before reclamation. 

Mill Run-B was opened in 1973. Prior to that the satellite data identified sev- 
eral pixels in the 1972 image in are;'.s where roads had been cleared in prepara- 
tion for strip mining. The 1973 aircraft and LANDSAT images indicate approxi- 
mately 65 acres of land had been affected by stripping operations. The mine was 
backfilled and seeded during the next year, reducing the 1974 acreage to 53 acres. 
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The mines on Franklin Hfll (Figure 2) were operative throughout the period 
covered by this Investigation. iUgh altitude photography ai the region In 1972 
and 1974 shows that these mines did not alter slgiiiflcantly in surface area over 
that period, even though 100 thousand tons of coal were extracted from these 
mines in 1973 alone. The Fronklin Hill mines were opened befcie iS72 and 
the same areas have been reotnpped several times. The Increased area Iden- 
tified in the 1974 LANDSAT image is due to new spoil piles and enlarged haul 
roads. 

Figure 9 shows the band-ratio classifications applied to the same study area 
as Figure 6. The results demonstrate the geographic cxtendibility of the ratio 
signature. False-alarms are similar to those resulting from the 4-channel 
classificatitja. ^Q8classificatlon of this sort are not generally serious since 
the locationh of most of the mines are well known. 

CONCLUSIONS 


This sbidy has demonstrated the feasibility of strip mine monitoring with 
LANDSAT multlspectral data to within 2 hectares (5 acies). The average 
acc'.racy ol clasriflcatlon is greater than 939f for LANDSAT images from three 
dai'. .. Using band-ratioing techniques it is possible to extend signatures over 
a large g>iographic area andtcmporall* *^o other LANDSAT images. 

The procedures developed in this study could be incorporated into a compre- 
hensive monitoring program to provide, in a rapid and inexpensive manner, 
accurate information on the location, size, and condition of areas affected by 
surface mining. Multlspectral analysis of satellite digital data would be uscfid 
in validating opcip.tors’ reports on the size and status of mining operations, 
locating abandoned and unrevegetated mines, aad assessing reclamation costs 
and requirements. 
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ngure 1. LANDSAT Mosaic of Maryland with the Test Site and I^arger Study Region 
Outlined; Insert Map of Coal Basins (from Yokes and Edwards, 1968). 



Figure 2. Low Altitude NASA/Wallops Aircraft Photo-Mosaic, Taken iO/lG/73, 
Showing the 51.8 Sq. Km Test Site with Test Strip Mines Identified. 


14 



Figure 3. Computer Output from Video Monitor Showing Test Site 
from 1 September 1973 Scene (1405-15242). 


15 



AIRCRAFT MEASUREMENTS; HECTARES 



LANDSAT MEASUREMENT; HECTARES 

(1 HECTARE - 2.47 ACRES) 


Figure 4. Linear Regression Analysis of MSS 5/6 Strip Mine Classificatioii. 
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Figure 6. Four-Band Classification of Strip Mines and Partially 
Revegctated Mines in Entire Study Region. 






Figure 7. Test Site Strip Mines Identified on Three LANDSAT Images from 
Same Spectral Signature in Ban<l 5/6 Ratio. 



i 


6 SEPT 1972 




22 JULY 1974 




Hgure 8. Strip Mines Identified on Three LANDSAT Images from Same 
Speetral Signature In Band 5/6 Batio. 
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Figure 9. Ikind- Ratio C'lassification of Strip Mines for the Study Region, 

1 September 1973. 


